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We study massive particle radiation from cosmic string kinks, and its observability in extremely 
high energy neutrinos. In particular, we consider the emission of moduli — weakly coupled scalar 
particles predicted in supersymmetric theories — from the kinks of cosmic string loops. Since kinks 
move at the speed of light on strings, moduli are emitted with large Lorentz factors, and eventually 
decay into many pions and neutrinos via hadronic cascades. The produced neutrino flux has energy 
E > 10 11 GeV, and is affected by oscillations and absorption (resonant and non-resonant). It is 
observable at upcoming neutrino telescopes such as JEM-EUSO, and the radio telescopes LOFAR 
and SKA, for a range of values of the string tension, and of the mass and coupling constant of the 
moduli. 
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I. INTRODUCTION 

Theories with spontaneous symmetry breaking usu- 
ally have topologically non-trivial vacuum configurations. 
Depending on the topology of the vacuum after the sym- 
metry breaking, stable relics called topological defects 
— such as monopoles, strings or domain walls — could 
be formed in the early universe [T]. Strings can form if 
the vacuum manifold is not simply connected. Although 
monopoles and domain walls are generally problematic 
for cosmology, cosmic strings are compatible with the 
observed universe, provided that their tension is not too 
large (Sec. [n] see, e.g., Refs. [2HD] for reviews). Cosmic 
strings are predicted in grand unified theories (GUTs) 
and superstring theory, and their existence can be re- 
vealed through their effects on the cosmic microwave 
background (CMB), large scale structure and 21 cm line 
observations, and — more directly — by detecting their 
radiation, such as gravitational waves and cosmic rays. 

Since cosmic strings have GUT or superstring scale en- 
ergy densities in their core, they can be significant sources 
of ultra high energy {E > 10 11 GeV) cosmic rays [7J- 
113] . either as isolated objects, or possibly in combination 
with other topological defects, like in monopole-string 
bound states [TlrfTB] . Among the cosmic rays, neutrinos 
are especially interesting. Their weak coupling to matter 
makes them extremely penetrating, so they are the only 
form of radiation (together with gravitational waves) that 
can reach us from very early cosmological times, namely, 
all the way from redshift z ~ 200 (see Sec. IV B). More- 
over, in the spectral region of interest, E 10 11 GeV, 
the neutrino sky is very quiet, since this region is be- 
yond the range of neutrinos from even the most extreme 
hadron accelerators (gamma ray bursts, supernova rem- 
nants, active galactic nuclei, etc.). Therefore, even a 
low statistics neutrino signal beyond this energy would 
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constitute a clean indication of a fundamentally different 
mechanism at play, such as a top-down scenario involv- 
ing strings or other topological defects. Experimentally, 
the technologies to detect ultra high energy neutrinos are 
mature: they look for radio or acoustic signals produced 
by the neutrinos as they propagate in air, water/ice, or 
rock. After the successful experiences of ANITA [19j . 
FORTE [20], RICE [3T] and NuMoon [22] — the lat- 
ter using radio waves from the lunar regolith via the so 
called Askaryan effect [25J — a new generation of exper- 
iments is being planned, that can probe neutrinos from 
cosmic strings with unprecedented sensitivity. Of these, 
the space based fluorescent light telescope JEM-EUSO 
[21] . and radio telescopes LOFAR [25^ and SKA (2D] seem 
especially promising. 

One of the distinguishing effects of cosmic strings as 
cosmic ray emitters is that they can produce bursts from 
localized features called cusps and kinks (Sec.|ll]), where 
ultrarelativistic velocities are reached. The radiation 
from cusps and kinks is very efficient, whereas the emis- 
sion from cusp/kink-free string segments is exponentially 
suppressed. This enhanced emission has been studied in 
connection with gravitational waves [2~TH2"9] , and electro- 
magnetic radiation [371 HH1 ESI EI] like gamma ray bursts 
[52H54] and radio transients [55H57J, as well as neutrino 
bursts [ID] . 

Among the several scenarios considered, there are a 
few that predict cosmic ray and neutrino fluxes at an ob- 
servable level, e.g., Refs. [IDHI5] . One of these, Ref. [15] . 
involves the decay of moduli — massive scalar fields that 
arise in supersymmetric and superstring theories — , that 
can have various masses and couplings to matter. Mod- 
uli with coupling stronger than gravity are fairly natural 
[55H12] and relatively unconstrained due to their very 
short lifetimes |43| . compared to gravitationally coupled 
ones [44 46 . By decaying into hadrons, the moduli even- 
tually generate a neutrino flux. In Ref. [13] the emission 
of such moduli from string cusps, and the corresponding 
neutrino flux were discussed. 

In this paper, we elaborate on the theme of moduli- 
mediated neutrino production from strings, and study 
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modulus emission from kinks. We show that the emission 
from kinks is very efficient, and is the dominant energy 
loss mechanism for the cosmic string loops for a wide 
range of the parameters. We calculate the neutrino flux 
expected at Earth after a number of propagation effects, 
mainly absorption due to resonant (Z° resonance chan- 
nel) and non resonant neutrino-neutrino scattering. We 
find that the flux might be observable at near future sur- 
veys, JEM-EUSO, LOFAR and SKA, depending on the 
parameters. 

The structure of the paper is as follows. After dis- 
cussing some generalities on strings and kinks in Sec. [XXJ 
the modulus emission from a cosmic string kink is cal- 
culated in Sec. |ffl| In Sec. llVj we discuss the decay of 



moduli, the properties of the hadronic cascade initiated 
by their decay into gluons, and propagation of extremely 
high energy neutrinos in the universe. In Sec. [Vj es- 
timates are given for the kink event rate, the neutrino 
flux, and its detectability by the existing and future neu- 
trino detectors. We also discuss the constraint from high 
energy gamma ray observations. Finally, in Sec. VI we 
give our conclusions. 



II. COSMIC STRINGS 

Much of the phenomenology of a cosmic string depends 
on its tension (or mass per unit length), \x. It is often ex- 
pressed in Planck units, as G/i, where G is the Newton's 
constant. Several cosmological and astrophysical obser- 
vations place upper limits on G/i; we briefly review them 
here. 

Since cosmic strings can create scalar, vector and ten- 
sor perturbations, they were initially considered as seeds 
for structure formation [47U49] , Hence, they contribute 
to the anisotropy of the cosmic microwave background 
(CMB) [SUM] and B-mode polarization of the CMB [581- 
[r?Tj . However, the current measurements of the CMB 
anisotropies at small angular scales by WMAP [62] and 
SPT [63], reveal that the cosmic string contribution to 
the total power is less than 1.75%, which translates into 
a constraint on the string tension, G/i < 1.7 x I0 -7 [57] , 

Although their contribution to the density perturba- 
tions is small, strings can still have effects on the early 
structure formation [64] . early reionization due to early 
structure formation [651 166] . formation of dark matter 
clumps [5Zj, and might yield detectable signal in the 21 
cm measurements [68H72] . Cosmic strings also produce 
gravitational waves 73 in a wide range of frequencies, 
both as localized bursts and as stochastic background, 
which can be detected by LIGO, eLISA and pulsar tim- 
ing array projects [29 ] l74Tf77] . The most stringent bound 
comes from the pulsar timing measurements, which put 
an upper bound on the long wavelength stochastic gravi- 
tational wave background, h 2 ficw ^ 5.6 x 10~ 9 yielding 
the constraint G/i <4x 10 -9 [74]. However, this upper 
bound is obtained by ignoring the kinetic energy of the 
cosmic string loops, and by assuming that cosmic strings 



only decay by emitting gravitational waves. Thus, the 
pulsar timing bound is expected to be somewhat relaxed 
by taking these effects into account. 

Cosmic string loops can emit moduli efficiently in the 
early universe when the length of the loop is of the order 
of the Compton wavelength of the emitted particle [44j . 
If moduli are gravitationally coupled to cosmic strings, 
very stringent cosmological constraints can be put on the 
string tension, G/i, and the mass of the modulus, m [4~4T 
[46] . On the other hand, if their coupling is stronger than 
gravitational strength, modulus radiation becomes the 
dominant energy loss mechanism for the loops, and the 
lifetime of moduli becomes a lot shorter. These relax the 
cosmological constraints on moduli significantly [43]. In 
this paper, we shall adopt the parameter space consistent 
with all the constraints mentioned above. 

Cosmic strings are born as smooth objects, but after- 
wards they undergo crossings and self crossings, which 
lead to truncations and successive reconnections. Every 
crossing produces a kink on the string after reconnection. 
The result of such processes are string loops with a few 
kinks [78] . Kinks are discontinuities in the vector tan- 
gent to the worldsheet characterizing the string motion, 
and gravitational and particle radiation is very efficient 
at kinks yielding waveforms with power law behavior in 
the momenta of the emitted particles [2"Tl |2"9"] . There 
are also transient features on loops called cusps, where a 
part of the string doubles on itself, that reach the speed 
of light momentarily. Cusps also produce radiation in 
bursts, with waveforms that have a similar power law be- 
havior. On the other hand, radiation from cosmic string 
loops with no cusps or kinks is exponentially suppressed, 
leaving the kink and cusp radiation as an interesting win- 
dow on the observable effects from strings. In the next 
section, we shall study massive particle radiation from 
cosmic string kinks. 



III. MASSIVE PARTICLE RADIATION FROM 
KINKS 

The free part of the action has Nambu-Goto term for 
a string of tension /x, and the massive scalar field term 
for the modulus of mass m 
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where g is the determinant of the spacetime metric g^ v 
and 7 is the determinant of the induced metric on the 
worldsheet, V(cr,r), given by j ab = g^X^X^. The 
interaction Lagrangian for the modulus field and the 
string has the form 
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where a is the modulus coupling constant, \x is the string 
tension and m p is the Planck mass. Ignoring back reac- 
tion effects, the equation of motion for the worldsheet in 
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the flat background g^ u = rj^u = diag(— 1, 1, 1, 1), and in 
the conformal gauge, where cr° = r and a 1 — a, is 



X-X" 



0, 



(3) 



with the gauge conditions X • X' = and X 2 + X' 2 = 1. 
The general solution can be obtained in terms of the right 
and left moving waves as 

X(a,r) = i[X + (a + )+X_(<7_)], (4) 

where a± = a ± r, and the gauge conditions are now 
given by 



X' ± 2 = 1, 



(5) 



where the prime refers to the derivative with respect to 
the corresponding light cone coordinate a±. 
The total power of particle radiation is 



P 



P, 
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where the power spectrum P n can be calculated by using 



dP„ 
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where G — m~ 2 is the Newton's constant, a is the modu- 
lus coupling constant, k is the momentum of the emitted 
particle, ui n — Airn/L = \Jk 2 + m 2 is the energy, L is 
the loop length, and T(k,u> n ) is the Fourier transform of 
the trace of the energy momentum tensor of the cosmic 
string loop given by 
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Using Eq. Q and the lightcone coordinates a±, Eq. Q 
can be factorized as 



derivative of the phase vanishes — or kinks — points 
where the vector tangent to the worldsheet has a dis- 
continuity. The cusp case has been studied for massive 
particle emission in Ref. [13]. At the cusp both $± are 
saddle points, hence their derivatives with respect to the 
corresponding lightcone coordinates vanish. On the other 
hand, for the case of a kink, either <E> + or has a saddle 
point, and the other one has a discontinuity. In what fol- 
lows we assume that $+ has a saddle point and $_ has 
a discontinuity. Then, assuming the kink is at a± = 0, 
worldsheet can be expanded about a± = as follows: 
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Using the gauge conditions one can show that 
&i, fi2 and X^ 1 ' are unit vectors, Xj/ • X^ 2 "* = and 
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The curvature of the string can 



be approximated as |X\ | ~ 2tt / L if the string is not 



too wiggly. Using the expansions (10) and (11), and the 
gauge conditions ([5]), the phases <E>+ and <&_ can be ob- 
tained as 
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for d— > 0, (13) 



where s\, S2 are constants of order 1, |k| = k and we 
assumed that k // X^\ It can be shown that [2"5h3T] 
when moduli are emitted at a small angle rather than 
being parallel to the direction of X + at the saddle point, 
the expansion (10) still applies provided that the angle 
satisfies 
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T(k,w„) = [ L da+ f L da-{l+X' + • X'Je* 

(9) 

where <I>± = cj„ct± =F k • X±. 

The integral in Eq. ^ is exponentially suppressed 1 
for a smooth loop of cosmic string of length L >> 1/m 
|44j . However, the phase become stationary if the string 
has cusps — saddle points on the worldsheet where the 



otherwise, leading to exponentially suppressed power. 
The term in the integrand of Eq. (|9| can be found as 
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where c and d are constants of order 1, which we will 



take as 1 in what follows. Using Eq. (15), to the leading 
order we obtain 



T(k,uj n ) 
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1 L 1/m condition is satisfied in the very early universe when the 
cosmic string loops can be quite small. Hence, moduli produced 
efficiently at earlier epochs can have cosmological effects, and 
subject to constraints I43H46| 
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These integrals can be written explicitly by using 
Eqs. (12 1 and (13) as follows 



da+e 
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After a change of variables, one obtains [T3] 

i/2 r°° ;3 ( , x 3 \ 
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where u = Lk — l) 3 ^". The imaginary part of the 
integral vanishes, and the real part is given in terms of 
the modified Bessel function of order 1 /3 
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The function if 1/3 (u) exponentially dies out at large u, 
and it can be approximated as a power law in the limit 
u « 1 as K\/^(u) ss u -1 / 3 . This limit corresponds to 
k >> m, and in this regime, we can write u ~ Lm 3 /16k 2 . 
Then, we obtain 
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where this formula is valid when u < 1, i.e., k > k c , 
where 
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For smaller values of k, I + is exponentially suppressed, 
thus we are only interested in the above regime for prac- 
tical purposes. 



Using Eqs. (13) and (17), the integral i_ can be simi- 
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which results in 

/_ ~ ^k~\ 
where the sharpness of a kink is defined as 

if) ee -(1 - fii • n 2 ). 



(23) 
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Using Eqs. \2l\ and ( |24[ ), we find the power spectrum 
from Eq. Q as 

ePP iba 2 Gn 2 , /n, 9/q , ^ 1 / — - . . 

L ' k ' , k>-mVmL. (26) 
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Integrating over solid angle gives a factor 

Q k - 27r6» fc - 2n(kL)- 1/3 , 



where 9 k given by Eq. (14) is used. Then, the total power 
can be obtained as 
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where we used the cutoff for upper limit for the integral 
over momenta [T2l l80l 



3/4^1/2 
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and the lower limit fc m i n ~ k c from Eq. (22). Note that 
for typical values of the modulus mass m and the string 
tension fi, the logarithmic factor is about 20. Then, we 
can simply write the total power as 
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where we define a ee y/ipa. Number of particles emitted 
from a kink with momenta k in the interval (fc, k + dk) 
can be found from Eq. (p6| as 
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In addition to moduli, cosmic string loops also produce 
gravitational radiation with the power [2] 

P g - 50G/A (32) 



It is convenient to define the power as 
P = TG/i 2 , 



where 



a 
■50 



for a 2 > 50 
for a 2 < 50. 



(33) 
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The dominant energy loss mechanism for loops deter- 
mines the lifetime of a loop as 



tl 
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Then, the minimum loop size that survives at cosmic 
time t is 



TGfit. 



(36) 



In the next section, we shall discuss the decay of the 
moduli produced from cosmic string kinks into neutri- 
nos via hadronic cascades, and the propagation of these 
neutrinos in the universe. 



IV. PARTICLE DECAY AND PROPAGATION 

For simplicity, throughout this paper we will assume a 
matter dominated flat universe model, which lets us carry 
out the calculations analytically. We assume cosmolog- 
ical constant A = 0, and the total density parameter 
O m + fi r = 1 has matter and radiation components. We 
use the following values of the cosmological parameters: 
age of the universe to = 4.4 x 10 17 s, time of radiation- 



(27) matter equality t cq = 2.4 x 10 12 s, 1 + z c 



3200 



The scale factor in the radiation and matter eras are 
respectively given by a T oc i 1 / 2 and a m cx t 2 / 3 . Using 
a /&o = 1/(1 + z )i a o = 0) t ne cosmic time can be written 
in terms of redshift as t = ta{l + z oq ) 1 / 2 (l + z)~ 2 in the 
radiation era, and t = to(l + z)~ 3 l 2 in the matter era. 
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A. Modulus decay 

The decay channel for moduli with the largest branch- 
ing ratio is the decay into gauge bosons with the interac- 
tion of the form 139 



C 



(37) 



for a modulus field <f> and a gauge field of field strength 
Fa V - For the gauge bosons in the standard model, the 



modulus lifetime is estimated as 



t ~ 8 x 10" 6 to; 



s. 



(38) 



Since most of the moduli are emitted from kinks with 
momenta k ~ my/mL/A [because of the decreasing power 



law given by Eq. ( 31 1] , their lifetime is boosted by a factor 
of 7 ~ s/mL/i. For the fiducial values of the parameters, 
the Lorentz factor of a modulus emitted at redshift z and 
survive at present epoch is given by 



7(2) 



^mTGfit 
4(1 + z) 



6.5 x 10 
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where we have used the fact that loops of size L m [ n given 



by Eq. ( 36 ) yield the dominant contribution to the ob- 
servable events, the factor of (1 + z) in the denominator 
takes into account the redshifting of the energy of the 
moduli emitted at epoch z. Thus, the ratio of the lifetime 
of a modulus emitted at redshift z < z oq and decaying at 
redshift z d to the cosmic time at epoch z^ is 



t{z,z a ) 
t(z d ) 



10 



-3 



1/2 
M-17 
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r i/2 m 5/2 (1 + 3)7/4 
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Note that T ^ 50 from Eq. (34) and z A ^ z. Hence, 

as they 



< 



moduli will decay in the same epoch, z d 
are produced. Therefore, we assume that all the moduli 
decay before they reach the Earth. 

The most efficient channel for neutrino production 
from modulus decays is the decay into gauge bosons. In 
particular, gluons decaying into hadrons produce neutri- 
nos with the largest multiplicity [TU1 13] • The interaction 
of a modulus with a gluon field is of the form (37), and 



the hadronic cascade from these gluons produces numer- 
ous pions of either sign, which eventually decay into neu- 
trinos and antineutrinos. For both, we expect a flavor 
composition in the ratio u^:v e :v T = 2:1:0, from the 
pion decay chain. 

The number of neutrinos per unit energy can be found 
by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) method. Monte Carlo simulations for the 
hadronic decay of a very massive particle show a power 
law behavior in energy as E~ n with n = 1.9 [SI]. For 
simplicity, we approximate the index as n ~ 2. Then, 
the fragmentation function has the form [10] [13] 



dN v 
dE 
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(l + z)E 2 
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where fc/(l + z) and E are the modulus and neutrino 
energies at the present epoch respectively. Here E m i n < 



E < E max jT3j, where 
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e-y(z) - 6.5 x 10 



13 



pl/2, 1/2 1/2 
(1 + z) 7 / 4 



GeV, 
(42) 

and E max ~ 0.1 k. We take e GcV = e/(lGeV) ~ 1 [13] . 
Since the neutrino spectrum has the form E~ 2 , most of 
the neutrinos will have the energy E ~ £min- This in- 
troduces a lower bound on the redshift, below which no 
neutrinos are produced with a given energy E < .Emm- 
For our estimates, we are interested in energies E > 10 11 
GeV corresponding to the minimum redshift in the mat- 
ter era 



,(£,J~4U1 /i_ 17 m 4 n, n , 



(43) 



where En = E/fiO 11 GeV). Since the maximum red- 
shift from which neutrinos can propagate to us is set by 
the neutrino horizon z v <~ 200 (see Sec. IV B), requiring 



< 



z v , we have the constraint on the parameters 



Gn < 2.8 x 10" 15 r 



E 



ii- 



(44) 



B. Neutrino propagation 



The neutrino flux at Earth is affected by a number of 
propagation effects: the redshift of energy, flavor oscilla- 
tion, quantum decoherence and absorption. The redshift 
of energy will be included as we carry out the flux cal- 
culation in the next sections; the other effects, instead, 
warrant a separate discussion, which is the subject of this 
section. 

The oscillations of very high energy neutrinos have 
been discussed in detail (see, e.g., [H2]). Oscillations in 
vacuum are a good approximation, as the refraction po- 
tentials due to the intergalactic gas and to the cosmolog- 
ical relic neutrino background (which is assumed to be 
CP-symmetric here) are negligible [S3J. For the large 
propagation distances we consider, the flavor conver- 
sion probabilities are energy independent, as the energy- 
dependent oscillatory terms average out [84 . For our 
predicted initial flavor composition, Vu:v e :v T — 2:1:0 



(Sec. IV A), the effect of oscillations is to equilibrate the 



flavors 82 , therefore the composition at Earth should be 
v^:v e :v T = 1:1:1, for both neutrinos and antineutrinos. 

A neutrino oscillates as long as its wavepacket remains 
a coherent superposition of mass eigenstates. Depend- 
ing on the size of the produced wavepacket, decoherence 
can occur as the neutrino propagates, due to the differ- 
ent propagation velocities of the mass states. Dedicated 
analyses 84, 85] have shown that neutrinos of the ener- 
gies of interest here remain coherent over cosmological 
distances, therefore we do not consider decoherence ef- 
fects. 

Absorption effects are largely dominated by scattering 
on the relic cosmological background [SHI [57] , with neg- 
ligible contribution from other background species. In 
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first approximation, absorption can be modeled as a sim- 
ple disappearance of the neutrino flux; secondary neutri- 
nos generated by scattering are degraded in energy and 
therefore they are negligible compared to primary flux. 

The survival probability for the primary neutrinos, of 
observed energy E and production redshift z, is defined 
as [86H89]: 



P{E,z) = e- r "< B '*> , (45) 
where the optical depth for the relic neutrino background 



is 



and 



t v (E,z)= dt' a vv {E,z)n v (z), (46) 

Jt(z) 



dt 1 



dz 



2 (1 + z) 5 / 2 ' 



(47) 



in the matter era. Here n v (z) = 56(1 + z) 3 cm~ 3 is 
the number density of relic neutrinos in each of the six 
species (neutrinos and antineutrinos of each flavor), and 
a vv (E,z) is the neutrino- neutrino cross section, eval- 
uated at the production energy E' — E(l + z), and 
summed over all the neutrino species in the background. 
For the energies of interest here, and at the leading order, 
this cross section is the same for neutrinos and antineu- 
trinos, and is practically flavor-independent: 

= er e — ct m ~ cr T (48) 
a a = 2J + -> any) + o(v a + Dp -> any)} . 

In the limit of massless neutrinos, m v m 0, the Z°- 
resonance effects can be ignored and the maximum cross 
section is attained at E > 10 11 GeV [13] : 



N 



TT 



(49) 



where TV - 10 - 15, G F = 1.17 x 1(T 5 GeV 2 , and 
m w ~ 80.39 GeV. Using Eq. Q in Eq. ((46}, and re- 



quiring r„ — 1 for absorption, the neutrino horizon — the 
maximum redshift from which the neutrinos with energy 
E can propagate to us — is given by [SHI HI] 



200, 



(50) 



for energies E > 10 11 GeV. In this regime, P(E, z) can 
be approximated as a step function 



P{E,z) « 1 -9(z-z„), 



(51) 



which becomes handy when estimating the neutrino flux 
analytically. 

If the Z° resonance is realized in the annihilation chan- 
nel v a + v a —¥ any, at a redshift z rcs along the neutrino 
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FIG. 1: Gr aph ical representation of the interval of observed 
energy, Eq. ( 54 1, where suppression of the flux due to resonant 
v — v annihilation is expected (region between the diagonal 



lines) . The dashing indicates the region where Eq. ( |54[ ) is only 
indicative, due to thermal effects influencing the resonance. 
The horizontal shaded area refers to the interval of neutrino 
masses where the neutrino mass spectrum is strongly degen- 
erate. We also mark the value m„ ~ 0.05 eV, which is the 
highest mass expected for non-degenerate (hierarchical) mass 
spectrum. 



path, a pronounced dip in the neutrino spectrum is ex- 
pected at the resonance energy due to the strong en- 
hancement of the cross section E 



{s/Ml 



1) 



r|/M 



(52) 



where s ~ 2m u E(l+z rcs ) if the background neutrinos are 
not relativistic. The effect of the resonance is especially 
transparent in this case [87l |88] ; we discuss it here in its 
essentials. 

Considering their momentum, p v {z) ~ 6.104x 10~ 4 (1+ 
z) eV, the cosmological neutrinos are non-relativistic to- 
day for masses exceeding ~ 10~ 3 eV, and throughout the 
interval of redshift of interest, z ^ z v , if nij 3> p v (Zu) — 
0.1 eV. From the data of oscillation experiments (see, 
e.g., }92j ) it is known that, above this value, the neutrino 
mass spectrum becomes degenerate: mi ~ m 2 — m 3 . 
Therefore, we can reason in terms of a single neutrino 
mass value, m v , and take m v = 0.3 eV as reference. The 
degenerate case is optimal for the observability of the res- 
onance effect, because the dip in the spectrum occurs at 
the same energy for all neutrinos and has a sharp shape. 
Furthermore, it is located in the region of the spectrum, 
~ 10 — 10 13 GeV, where experiments have good sensi- 
tivity [Ml (see Fig.[l). 

For a neutrino of energy E at Earth, the Z° resonance 
is realized at redshift z res if 



E = 



1m v (1 + Z res ) 



1.4xl0 13 GcV 



0.3 eV 



(1 + Zrcs) 1 
(53) 
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with M z ~ 91.19 GeV the mass of the Z° boson. It 
follows that the flux of neutrinos of observed energy 

E = (6.9 x 10 10 - 1.4 x 10 13 ) GeV , (54) 

is affected by the resonance between z = z v and the 
present epoch (see Fig. [I]), and therefore should be 
strongly suppressed compared to the flux at energies out- 
side this interval, where the smaller, non-resonant, ab- 
sorption cross section is at play. 

Following the detailed discussion in Ref. [55], we cal- 
culated P(z,E) and used it to obtain the neutrino flux 
expected at Earth from all sources at all redshifts. This 
flux is calculated by convolving the flux per unit of pro- 
duction redshift with the probability P(z, E); it exhibits 
the characteristic suppression dip in the interval given in 
Eq. (54), as expected (see Sec. VC). 



The absorption pattern is more complicated if the neu- 
trino mass spectrum is not degenerate, i.e., uii ^ m 2 <C 
TO3 ~ 0.05 eV (or, m 3 <C mi ^ m 2 — 0.05 eV). For this 
configuration the probability P(z, E) has three distinct 
dips of resonant suppressions at separate resonance en- 
ergies [531 [55], corresponding to the three masses. These 
dips are broadened, in energy, by the integration over the 
production redshift, and, most importantly, by thermal 
effects, which are important for in this range of neutrino 
masses [HH [55J [5T] . We postpone a discussion of these 
effects to a forthcoming publication 93 . 



V. NEUTRINO FLUX AND DETECTION 

As kinks move along a loop of cosmic string, they emit 
particles in a fan-like pattern, and scan a ribbon of solid 
angle fl ~ 2^6*^ [see Eq. (14)]. Thus, one can analogously 
visualize the radiation from a kink as a source of light 
emitted from a lighthouse passing by. An observer who 
happens to be within the beam direction sees particles 
as a burst event provided that the flux is detectable. In 
this section, we make order of magnitude estimates for 
the event rate for bursts and neutrino flux, and compare 
it with the existing and future neutrino experiments. 



A. Loop distribution 

The distribution of cosmic string loops has been stud- 
ied both analytically [51H55] and in simulations [55HTU7] . 
Although there seems to be a consensus on the distri- 
bution of subhorizon size large loops, there is still not 
a good understanding of the small loop distribution. In 
what follows, we use the results from the latest simula- 
tion that has the largest dynamical range up to date for 



the evolution of the cosmic string network |107j . where 
it has been confirmed that the large loops form with size 
(3t, where /3 ~ 0.1, and t is the cosmic time at which the 
loop is chopped off the network of long cosmic strings. 



The density of long strings is p ~ Ca*A 2 > with ( ~ 16. 
Using this framework, we can estimate the number den- 
sity of loops of length (L,L + dL) that are formed in the 
radiation era and still survive in the matter era as 

n(L, t)dL ~ p- 1 C( ( S* e q) 1/2 t" 2 L- 5/2 di, (55) 

where TGi-d L 0t eq and p is the reconnection prob- 
ability. There are also loops formed in the matter era, 
however, we have verified that their number density is 
negligible compared to the loops surviving from the ra- 
diation era given by Eq. (55). The dependence of loop 



density on reconnection probability has not been resolved 
yet, however, it is expected that the loop density in- 
creases for decreasing reconnection probability as dis- 
cussed in Refs. [55J I108j . For ordinary cosmic strings, 
p = 1, and it has been estimated as 10~ 3 ^ p ^ 1 for 
cosmic F- and D-strings [109] , Note that the most nu- 
merous loops have size of order L min ~ TGpt. As we 
shall see in Sec. [Vj those will give the most dominant 
contribution to the observable effects, such as the diffuse 
neutrino flux. 

When a loop of size f3t eq is formed, it will decay by the 



time [see Eq. (35 )] 



fit 



eq 



(56) 



This loop can survive until epoch z < z cq provided that 



(57) 



i + z>o.07r 2/3 ^ 2/ 8 3 . 



where = G/x/10 -8 , T > 50 is given by Eq. |34| , 
and we used t = to (1 + z)~ 3 ^ 2 . Hence, we can conclude 
that even for the maximum G/i allowed by the current 
bounds, the loops can survive all the way to the very 
recent epochs from which we can get observable effects 
unless T is too large. 



B. Burst rate 

The number of kink bursts per unit time can be esti- 
mated as [T01IT5] 



, . dtt n(L,t)dL „ T , . 

dN = ; ■ ' \ , dV(z), 

4ttL 1 + z /2 v h 



(58) 



where n(L,t)dL/(L/2) is the frequency of a kink event 
per physical volume per unit loop length, L/2 is the oscil- 
lation period of a loop of length L, dVl/Aii ~ Ok/2 is the 
probability that an observer lies within the solid angle of 
kink radiation, and dV{z) is the physical volume in the 
interval (z, z + dz) in the matter era, given by 

dV{z) = 54^ 3 (1 + zr 11/2 [(l + z) 1 ' 2 - l] 2 dz. (59) 

To find the total burst rate, we integrate Eq. ( [58] ) over 
L and z. Integral over L is dominated by its lower limit 



L m { n given by Eq. ( 36 ) , and the integral over redshift is 
dominated by its upper limit z v ~ 200. Numerically, we 
obtain the total event rate as 



N ~ 1.6 x 10 18 p- 1 T- 3 fiZli m 4 1/2 



yr 



(60) 



Remember that r > 50 is given by Eq. (34). Since ex- 
periments run for a few years, the event rate should be 
at least ~ 1 per year to get observable events. Requiring 
N > lyr -1 yields the constraint on the parameters 



Gn< 1.2 x io- 11 ^- 173 r _1 ' 



^4 "^yr ^ 



(61) 



where N v 



N/(lyr). 
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C. Neutrino flux 

The diffuse neutrino flux is obtained using the flux 
from a single kink on a loop, and summing over all the 
loops in a volume constrained by the neutrino horizon 
z v ~ 200. It can be estimated by 



ME) 



dN 



£,(E,z)dN(k)P(E,z), (62) 



where 



r(z) = 3*o(l + z)- 1/2 [(l + z) 1 ' 2 - 1], (63) 



is the physical distance to the source at redshift z in the 
matter era, dN is the kink event rate defined by Eq. ( 58 ) 
and fifc is the solid angle into which moduli are emit- 



ted given by Eq. (27). Here £(E,z) is the fragmentation 
function given by Eq. (J4TJ) , dN(k) is the number of mod- 
uli emitted from a kink with momenta k in the interval 
(k, k + dk) given by Eq. ( |3l| ), and P(E, z) is the survival 
probability of neutrinos defined in Eq. (45). Putting ev- 
erything together, we obtain 



e 2 me) ~ omp-^^c ( - 



1/2 

'eq \ -2//-I \2 2j.-l/2 

1 a (Gfi) m t 



dL 



(1 + Z)V* J Lmin L5/2 7 fcmin k 



dk 



.(64) 



Note that the integral over k gives a logarithmic factor 
ln(k max /k min ) - lnGtz 1 / 2 /m) 3 / 2 - 20 from Eqs. Q and 



(29), and the integral over L is dominated by its lower 
bound L m j n given by Eq. ( |36[ ). The integral over z can 
be done numerically. However, it is useful to see the lim- 
iting case m v ~ 0, where we can ignore the Z-resonance 
effects, and carry out the integral over redshift analyti- 
cally. Using the approximate form of P(E, z) given by 



Eq. (51) in Eq. (64), the neutrino flux can be calculated 
as 



E 2 M E ) - 2 x 10 



!/ 2 -2 

r 3 / 2 



GeV 

cm 2 s sr 
(65) 



FIG. 2: Examples of neutrino flux from cosmic string kinks, 
via moduli decay. We took G/i = ICP 17 . The dashed curves, 
(a),(c) and (d) refer to m„ = and a — 1 , 10 2 , fO 3 respec- 
tively. The dot-dashed line, (b), is the same as (a) but for 
m v — 0.3 eV, with resonant absorption effects included. Cases 
(c) and (d) are not affected by the Z° resonance for this value 
of the neutrino mass. The figure also shows existing limits 
from ANITA, FORTE, NuMoon and RICE [T9h22]. and ex- 
pected sensitivities of the future detectors JEM-EUSO, LO- 
FAR and SKA I24H26]. 



Using z m i n from Eq. ( 43 1 , the predicted diffuse neutrino 
flux in the m v ss limit is 



E 2 M E ) ~ 8x10 



3/7 - 2 it,1/7 
-5 A*-17 a E ll 



P 



pll/7 



[l-(^min/^) 1/4 ] 



GeV 



cm' s sr 
(66) 

Taking into account the neutrino mass in the survival 
probability P(E,z), and fixing m„ = 0.3 eV and recon- 
nection probability p — 1, we evaluate Eq. (64) numer- 
ically. In Fig. [2] we show the predicted flux for a few 
different values of the parameters G\i and a, together 
with the detectability limits of the current and future 
neutrino experiments. 



D. Diffuse gamma ray background constraint 

As moduli decay via hadronic cascades, the pions from 
this process also decay into photons and electrons. These 
high energy photons and electrons interact with the CMB 
photons and extra galactic background light, producing 
an electromagnetic cascade, whose energy density is con- 
strained by the measurements of diffuse gamma ray back- 
ground |110j . The strongest upper bound on the cascade 
energy density comes from the highest energy end of the 
observed spectrum. The most recent data from Fcrmi- 
LAT observations reach E ~ 100 GeV |lllj . The cascade 
photons with energy E > -E a bs will be strongly absorbed 
due to interaction with the CMB photos, where E a ^ s due 
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to pair production can be estimated as 



-Eabs(z) 



£cmb(1 + z) 



5.6 x 10 5 



1 + 2 



GeV, (67) 



detector 2 . Time coincidence at arrival is expected for 
neutrinos from a single burst, because the emission oc- 
curs with a very short, practically vanishing, time scale 
of order El [Ml, 



where e C MB = 2.35 x 1CT 4 cV and m e = 0.511 MeV. This 
implies that a cascade photon of energy above E cas ~ 
100 (1 + z) GeV, is efficiently absorbed at redshift 



fcmin')' 2 



10- 42 s, 



(72) 



l + z c 



£abs(0) 
£cas(0) 



1/2 



for the fiducial values of the parameters, where fc m i n and 
7 are respectively given by Eqs. (22 1 and (39). The time 



70. (68) lag due to the spread in neutrino velocities is negligible 

for the Lorentz factors of interest here, 7 > 10 11 [see 



The electromagnetic energy density of the cascade 
from cosmic string kinks is jTOJ [T3] 



dt 



n(L,t)dLdP(k), (69) 



where /„■ is the fraction of energy transferred to pions in 
the hadronic cascade initiated by a modulus decay, 1/2 
is the fraction of energy transferred to electron-positrons 
and photons from the pion decays, dP{k) is the power 
emitted from a kink, given by Eq. ( |26[ ). Integrating over 
L, k and z [similar to the diffuse flux in Eq. (64)], and 
integrating over z up to z cas ~ 70, we have 



!/ 2 -2 tt 

k M i7« eV 

J cas - 1.2 X 10- 5 2=^- 



cm° 



(70) 



where r > 50 is given by Eq. (34). 



The maximum value of w cas allowed by Fermi-LAT dif- 
fuse gamma ray data is uj™ x ~ 5.8 x 10~ 7 eV/cm 3 |112j . 



Therefore, uj cas < oj™f x is satisfied for 



G/i<3x 10" 



(71) 



Note that ui cas > w™f x i s n °t strictly ruled out. This 
bound only constrains the observed highest energy dif- 
fuse gamma ray photons at E 7 ~ 100 GeV that are orig- 
inated at very large redshifts z < 70. The constraints on 
the energy density of cascade photons produced at red- 
shifts larger than z cas ~ 70 are much more weaker since 
they are more efficiently absorbed. Besides, the radia- 
tion from kinks is not homogenous, but confined to be in 
a narrow ribbon of width 2ir8k <C 1. Unless the cosmic 
magnetic fields are strong enough, the beamed electro- 
magnetic radiation from cosmic string kinks might not 
diffuse efficiently, hence, the constraint might be relaxed 
significantly. Nevertheless, the examples given in Fig. [2] 
respect the cascade upper bound given by Eq. (71 ). 



E. Neutrino bursts from individual kinks 



Eq. (42)] 



The fluence of neutrinos with energy above E, from a 
kink on a cosmic string loop, can be estimated as [TUJ [13] 



F u (> E) = J 



({E,z)EdN(k)P{E,z) 
^V 2 ' 



(73) 



Using Eqs. ([27j, (|31j), (|41j), (|45| and (|63j), and integrating 
over momenta yields 



F v {> E) ~ 4.3 x 10" 



-2 15/4^/2^-1 
* M-16 1 ^11 



[ A /TT^-l] 2 (i + z)9/ 4 

where we take the loop length to be L m i n given by 
Eq. ( |36[ ) since these loops are the most numerous, as 
was discussed in Sec. |V C[ hence it is more likely to get 
a burst from such loops. 

We can now estimate how many neutrinos might be 
detected at a detector of effective area 



A. 



dct 



M 

m N 



(75) 



where the nucleon mass, mjv ~ 1 GeV, M is the tar- 
get mass, and <j v jq ~ 10 -31 cm 2 is the neutrino-nucleon 
cross section. The reference cross section is from recent 
calculations at E ~ 10 12 GeV (see e.g., [IT5HIT5] ) and 
is a reasonable approximation for higher energies as well, 
due to the slow rise of <j v n at these energies (less than 
oc E 1 / 2 ). A typical value of the target mass for neutrino 
detection is M ~ 10 21 g, which applies to JEM-EUSO in 
its nadir mode at energy E ~ 10 11 GeV (23]. 

We model a best case scenario by choosing the closest 
distance to the source, z ~ z m - m (neutrinos can only come 
from z m i n < z < z v ), and the regime a 2 > 50 (where 
r ~ a 2 ), for which the the number of emitted neutrinos 
is larger. 

The number of events in a detector due to a burst can 
be estimated as: 



K = F V {>E)A, 



del 



1 7^, 3 - 14 ,, 2 - 82 TP' 

1.7 a 3 v_ 17 E- 



0.86 

14 1 



^16 



(76) 



Before closing, we comment briefly on the possibility 
to identify the neutrino emission of individual kinks, i.e., 
bursts, rather than the diffuse flux. The signature of a 
burst would be two or more time-coincident events in a 



2 Due to the opening angle at emission, the coincident events would 
most likely appear in different regions of the detector field of 
view; this will allow to distinguish them from one another. 
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where A 16 = A det /(10 16 cm 2 ). The fact that N v > 1 
means that, for our parameters of reference, the identifi- 
cation of a burst by time coincidence of multiple events 
is possible in principle, although in practice instrumental 
backgrounds might be an obstacle. 

Requiring M v 1, implies a minimum value of G/i: 



Gfj, > 7.8 x 10" 



18^,-1.12^.-0.31 ,1-0.36 



a. 



14 



.4 



in 



(77) 



This has to be combined with the maximum value im- 

< Zv ~ 200 [see 



posed by the condition that z m - m 
Eq. |44|] : Gfi < 2.8 x 10~ 15 a^m^ 1 E\ A . The range of 



Gy, is further restricted by imposing that burst detections 
are frequent enough, say one per year at least. Using 

,-1/6^-1/3 



Eq. (61), this gives Gfi < 1.2 x 10" 17 a 



Therefore, the burst detectability is possible for a some- 
what narrow range of Gji. 

Note also that a detector's capability to see bursts de- 
pends on its energy sensitivity: for most of the parameter 
space, the neutrino emission is concentrated above the 
JEM-EUSO peak sensitivity, E ~ 10 11 GeV, therefore 
detection at JEM-EUSO might be hard. However, LO- 
FAR and SKA are expected to surpass the JEM-EUSO 
sensitivity at higher energies (see Fig. [5]), and therefore 
are more promising burst detectors. The detection of a 
burst would be an important signature of cosmic string 
kinks or cusps (see Refs. [TOl H3| for bursts from cusps), 
complementary to a possible diffuse flux observation. It 
would also help breaking the degeneracy between the two 
parameters, a and G/i, since the detected number of neu- 
trinos from a burst, (76), and the diffuse flux, (66), have 



different dependences on the parameters. Besides, even 
if only single neutrinos are detected, the rate of events, 
(60), can be used to help distinguish cosmic strings as 



the source, and break the degeneracy of the parameters. 



VI. SUMMARY AND DISCUSSION 

Cosmic strings loops form as a result of reconnection 
of long strings, and self-intersection of large loops. Kinks 
arise naturally as a result of these processes. We studied 
how kinks can radiate moduli, particles that arise in the 
supersymmctric models of particle physics, and that can 
have various masses and couplings to matter. The decay 
of moduli into pions via hadronic cascades produces a 
flux of neutrinos, which can be observable depending on 
the parameters. 

Specifically, we considered the string tension G/z, the 
modulus coupling constant a and mass m as free param- 
eters, and showed that neutrinos with energies E > 10 11 
GeV can be easily produced by cosmic string loops via 
this mechanism, with flux 



E 2 J V {E) ~ 1.7 x 10 



GeV 



3/7 - 2 r-,1/7 
-7 M-17 a -^11 

p(r/50) u /7 cm 2 ssr' 



(78) 



The hadronic cascade stops producing pions at the mod- 
ulus rest frame energy of order e ~ 1 GeV. In the rest 



frame of the loop, this energy is boosted by the Lorentz 
factor 7, so that the the minimum observed energy of the 
neutrinos is: 



E min (z) -4.6X10 1 



.(r/SOjVyffym^eGeV 

(1 + Z)V4 



GeV. (79) 



The neutrino flux is shown in Fig. [2] for representative 
sets of parameters; the termination of the flux at E m i n 
appears clearly. The figure also gives the flux sensitivity 
of various experiments, showing that the predicted flux 
is within reach for the next generation neutrino detectors 
such as JEM-EUSO, LOFAR and SKA. 

A distinctive feature of radiation from cosmic string 
kinks is that particles are emitted in a fan-like pattern, 
confined into a narrow ribbon, hence bursts from individ- 
ual kinks can possibly be identified by timing and direc- 
tional coincidence. In Eq. (76), we estimated the number 



of neutrinos emitted by a kink, and the corresponding 
number of events in a detector of a given effective area. 
We found that, for the fiducial values of the parameters 
used in our analysis, multiple neutrinos can be seen in 
the field of view of the detector. 

If ultra high energy neutrinos are observed at future 
experiments, what would be possible to learn? Top- 
down mechanisms would offer natural explanations, and, 
among those, cosmic strings would be a favored candi- 
date. Even in the framework of cosmic strings, how- 
ever, data analyses will necessarily be model-dependent, 
and various models would have to be considered. Our 
scenario involving moduli is a possibility among many, 
and other intermediate states leading to neutrino pro- 
duction are possible, e.g., modulus emission from string 
cusps [13 and heavy scalar particle emission from cusps 
of superconducting strings jlOj . Another possible gen- 
eration mechanism of extremely high energy neutrinos 
could be the KK mode emission from cusps and kinks 
of cosmic F- and D-strings. The emission of KK modes 
of gravitons from cusps was studied in Refs. |116l 1117] , 
and various cosmological constraints have been put on 
the cosmic superstring tension. Depending on the pa- 
rameters, observable neutrino fluxes might be produced 
by this mechanism as well. 

A discrimination between different models will require 
the combination of complementary data, probably the 
detection of gravitational wave/electromagnetic coun- 
terparts of neutrino signals |118[ 1119] , The identifica- 
tion of point-like sources of extremely energetic neutri- 
nos (bursts) would favor cosmic string kinks or cusps as 
sources, a hypothesis that would be substantiated fur- 
ther by the observation of accompanying gravitational 
wave and/or gamma ray bursts. To distinguish between 
kinks and cusps could be possible since the event rate is 
larger for kinks for the given values of the parameters. 

In addition to a possible discovery of topological de- 
fects, detecting a flux of ultra high energy neutrinos 
might reveal new pieces of the still incomplete puzzle 
of neutrino physics. Most interestingly, if the data show 
a Z° resonance dip, we might gather information on the 
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neutrino mass and have another, perhaps more direct, 
evidence of the existence of the cosmological relic neutri- 
nos. The information on the neutrino mass might be es- 
pecially important if at least one neutrino is light enough 
to evade a direct mass measurement in the laboratory. 

It is important to consider, however, that the extrac- 
tion of any information from data would be complicated 
by many theoretical uncertainties. Let us comment on 
the uncertainties and simplifying assumptions of our cal- 
culation. First of all, we worked in a flat matter dom- 
inated universe, and ignored the recent accelerated ex- 
pansion period of the universe, whose effect can be at 
most about a factor of a few in our final estimates. We 
also approximated the neutrino fragmentation function 
for the moduli decays as dN/dE oc E~ n , and used n = 2, 
whereas the numerical calculations yield n ss 1.9 |81j . 
In our estimates we take into account the reconnection 
probability p. For cosmic strings of superstring theory, 
namely, F- and D-strings |109j , p < < 1 , whereas for ordi- 
nary field theory strings p = 1. The flux, event rate and 
the chance of getting neutrino bursts is expected to be 
enhanced for cosmic superstrings with p < 1, compared 



to ordinary cosmic strings. We ignored the backreaction 
of modulus emission from kinks on the evolution of kinks. 
Since the total power from a kink is only logarithmically 
divergent [see Eq. ( 28 )], the effect of radiation is expected 
to smooth out the sharpness of a kink slowly. Finally, 
our treatment of the neutrino absorption due to reso- 
nant scattering on the neutrino backround is limited to 
relatively large masses, m v J; 0.1 eV, for which thermal 
effects on the background are negligible. The generaliza- 
tion to include these effects is forthcoming |93) . 
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